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Abstract 

Q-f We present a QCD analysis of the neutral current neutrino-nucleus interaction at the 

O small- x region using the color dipole formalism. This phenomenological approach 

is quite successful in describing experimental results in deep inelastic ep scatter- 
ing and charged current neutrino- nucleus interactions at high energies. We present 
theoretical predictions for the relevant structure functions and the corresponding 
j> ■ implications for the total NC neutrino cross section. 
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1 Introduction 
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The interaction of high energy neutrinos in hadron targets are an important 
probe to test Quantum Chromodynamics (QCD) and useful to understand the 
parton properties of hadron structure. The several combinations of neutrino 
and anti-neutrino scattering data can be used to determine the structure func- 
tions, which constrain the valence, sea and gluon parton distributions in the 
nucleons or nuclei. In a similar way as for charged-lepton deep inelastic scat- 
tering (DIS), the deep inelastic neutrino scattering is also used to investigate 
the structure of nucleons and nuclei. In the leading order quark-parton model, 
the structure function F 2 is the singlet distribution, F£ N oc xq s = xJ2(Q J rQ)i 
the sum of momentum densities of all interacting quarks constituents. These 
relations are further modified by higher-order QCD corrections. Currently, 
the theoretical description of experimental data in neutrino DIS is reasonable, 
where main theoretical uncertainties are the role played by nuclear shadow- 
ing in contrast with lepton-charged DIS, and a correct understanding of the 
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low Q 2 limit. However, nuclear effects are taken into account by using the 
nuclear ratios R = F 2 j AF 2 extracted from lepton-nucleus DIS, which could 
be different for the neutrino-nucleus case. The low-Q 2 region cannot be ad- 
dressed within the pQCD quark-parton model since a hard momentum scale 
Qo > 1 _ 2 GeV 2 is required in order to perform perturbative expansion. 

Concerning the important issues raised above, in Ref. [1] we have performed 
an analysis of small-x charged current (CC) neutrino-nucleus DIS using the 
color dipole formalism [2]. The structure functions F 2 N , xF£ N and the quan- 
tity AxF^ N were calculated and compared with the experimental data from 
CCFR [3-5] and NuTeV [6] by employing phenomenological parameterizations 
for the dipole cross section which successfully describe small- £ inclusive and 
diffractive ep DIS data. In addition, we predicted the nuclear ratios R 2 and R3 
and single out the size of nuclear effects in each case. It was found that small- a; 
data show geometric scaling property for the boson-hadron cross section and 
the structure function F 2 is in agreement with the phenomenological imple- 
mentation using the saturation models at the small- a; region [1] . The structure 
functions xF^ N was discussed in detail, where the sea content described by 
the quantity AxF£ N is well described. This investigation calls attention to the 
robustness of the color dipole formalism to describe the total neutrino cross 
section of ultra-high energy neutrinos at very high energies. 

Our goal in this work is to perform similar investigation for neutrino neutral 
current (NC) scattering, which is an ideal probe for new physics. As an addi- 
tional motivation to study NC neutrino interactions we quote the purpose of 
a new high-energy, ultra-high statistics neutrino scattering experiment, Nu- 
SOnG [7] (Neutrino Scattering on Glass). This experiment uses a Tevatron- 
based neutrino beam to obtain an order of magnitude higher statistics than 
presently available for the purely weak + e~ interactions and DIS events 
which is about two orders of magnitude larger than past measurements will 
also be obtained. The very high statistics will permit an electroweak mea- 
surement using the DIS data sample from NuSOnG, through the "Paschos 
Wolfenstein method" (PW) [8]. The best electroweak measurement using DIS 
events to date comes from the NuTeV experiment, which has observed an 
anomaly. The DIS measurement requires the knowledge of PDFs which de- 
scribe the momentum distribution of quarks as a function of Q 2 and can bring 
in theoretical uncertainties from sources such as nuclear effects and nuclear 
isospin violation. NuSOnG probably will improve statistics in comparison to 
previous highest statistics experiments NuTeV and CCFR. Theoretical uncer- 
tainties can be addressed by making a high statistics measurement of the PDFs 
on glass, in situ. NuSonG will generate an unprecedent sample of > 100M DIS 
events which can be used to measure six structure functions (three on neu- 
trinos and three for anti-neutrinos) as well as the strange and anti-strange 
parton distributions. The NuSOnG experiment intends to measure separately 
AF 2 = ^ if c (x, Q 2 )-F 2 NC (x, Q 2 ) in CC W ± exchange and NC 7/Z exchange 
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processes, which can help to constrain charge symmetry violation. 

In this work we present a determination of the small-a; structure functions for 
NC neutrino-nucleus within the color dipole formalism. This approach allows 
for a simple implementation of shadowing corrections [9] in neutrino-nuclei 
interactions. This paper is organized as follows. In Sec. 2, the structure func- 
tions F% N F£ N and the nuclear ratio Ra = F% A j ' AF^ N are investigated in 
this approach at small-x region, employing recent phenomenological parton 
saturation models. It is shown that at small-x, the NC boson-nucleon cross 
section should exhibit the geometric scaling property, that has important con- 
sequences for ultra-high energy neutrino phenomenology. Finally, we also an- 
alyze the charm contribution to the total NC neutrino-nucleus cross section. 
The prediction is consistent with the existing experimental determinations. In 
the last section we summarize the results and present our conclusions. 



2 Neutral current neutrino-nucleon interaction at high energies 

In the color dipole approach [2], the small- X deep inelastic scattering (DIS) 
is treated in terms of the interaction of the qq color dipole of size r with 
the target nucleon (or nucleus) which is characterized by the beam/flavor 
independent color dipole cross section. The color dipole formalism allows an 
all twist computation, in contrast with the usual leading twist approximation 
of the structure functions. In the neutral current neutrino interaction, the 
process Z° + N — > X is viewed as the result of the interaction of a color 
singlet dipole, in which the gauge boson Z° fluctuates into, with the nucleon 
target. The interaction is modeled via the dipole-target cross section, whereas 
the boson fluctuation in a color dipole is given by the corresponding light-cone 
wave function. The neutral current (NC) DIS structure functions [10-12] are 
related to the cross section for the scattering of transversely and longitudinally 
polarized Z° bosons. That is, 



where r denotes the transverse size of the color dipole, z is the longitudinal mo- 
mentum fraction carried by a quark and V'tl are the light-cone wavefunctions 
for the virtual Z° bosons with transverse or longitudinal polarizations. The 
small-x neutrino structure function F% N is computed from expressions above 
taking F 2 = F? + F^. Explicit expressions for the light cone wavefunctions 
in the charged current case can be found at Refs. [10, 11]. The corresponding 
expressions for the special case of massless quarks and considering an isoscalar 
target are presented in Ref. [12]. We quote Ref. [1] and Refs. [13, 14] for the re- 
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cent phenomenology using the color dipole picture in the case of CC neutrino 
interactions. 

For the Z° boson, the coupling to quarks contains both vector and axial- vector 
parts and it is convenient to consider the basis of helicity spinors to decompose 
the Z° light-cone wave function into distinct vector and axial- vector parts [15]: 

tfj 1 V, *, Q 2 ) = Vx lX2 (r, z, Q 2 ) - A x x ^(f, z, Q 2 ), (2) 

where Ai,A2 denote the quark helicities. The probability densities of quark- 
antiquark states |^f°| 2 and |^/f°| 2 are given by, 
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The vector and axial- vector parts of the Z° wave function are given by [15], 



= -^g f v ^2N c {±ie ±id r[z5 Xl , ± 6 X2 , T - (1 - z)6 Xl ^6 X2 , ± ]d r + ™A 1)± <W} 

A X ± \ X2 = V^9 f A^N c {-ie ±w r[z6 M , ± 6 X2tT + (1 - z)6 Xl ^5 X2 ,±\d r 

K (er) 
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A x ^ = ^g{^_l {5 Xl ,- X2 2\ 1 [2Q 2 z(l-z) + 2m)] - iS XlM e- i2X ^2m f d r } 

where, az = 4:iia em / sin 2 (26V) an d e 2 = z( 1 - 2;)Q 2 + mj. The axial-vector 
couplings are given by g{? = v /p(// — 2Qf sin 2 6V) and = ^fp(ll). Here, 
// and are the weak isospin and electromagnetic charge of quark /, re- 
spectively. The relative coupling strength of the neutral to charged current 
interactions is labeled by p, where p = 1 at tree level in SM. The weak 
mixing parameter, sin 2 9 W , is related at tree level to Gp, M z and a em by 
sin 2 6V = 47ra em / v / 2 G F M 2 Z . 

On the other hand, the dipole hadron cross section Odi V contains all informa- 
tion about the target and the strong interaction physics. Here, we consider two 
representative models for the dipole cross section. In what follows one takes 
the phenomenological parameterizations: (a) Golec-Biernat-Wusthoff model 
(GBW) [16] and (b) impact parameter Color Glass Condensate model (b- 
CGC) model [17]. Both models are able to describe experimental data on 
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inclusive and diffractive deep inelastic ep scattering at small-x. The b-CGC 
model is also successful in describing data on exclusive vector meson produc- 
tion at DESY-HERA and includes important aspects as impact parameter 
dependence, extended geometric scaling and saturation physics at very small- 
x. The idea behind it is to introduce the impact parameter dependence into 
the CGC model [18], where dipole cross section is given by: 



&di P (x, r)=2 J d 2 bN(x, r, b), 



N (x, r, b) 



3 -Aln 2 (BrQ sa t) 



rQ sat < 2 



(5) 
(6) 



where Y = ln(l/x) and Q sa ,t(x, b) is the impact parameter saturation scale 
defined as [17]: 



Q S3 .t(x,b) = ( — 
x 



exp 
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In the b-CGC model the evolution effects are included using an approximate 
solution to the Balitsky-Kovchegov equation. In particular, one has the slope 
Bcgc — 5.5 GeV~ 2 . The constants A and B are obtained from continuity 
conditions at rQ sat = 2. We quote the original papers for details on the pa- 
rameterizations and determination of their phenomenological parameters. Two 
important parameters are A and 7 S (saturation anomalous dimension), which 
determine the energy dependence of the saturation scale Q sat - Here, we use 
an effective light quark mass, to/ = 0.14 GeV and the charm (bottom) mass 
is set to be m c = 1.4 GeV (m b = 4.5 GeV). In the dipole cross section we use 
x = x Bj [l + (4m 2 f /Q 2 )]. 

For scattering on nuclei, we use the extension of the color approach for nuclear 
targets taking the Glauber-Gribov picture [9], without any new parameter. In 
this approach, the nuclear version is obtained replacing the dipole-nucleon 
cross section by the nuclear one, 



nucleus 

dip 



(x, r;A)=2jd 2 bl [ l- exp [-i T A (b) ^"(x, r)] } 



(8) 



where b is the impact parameter of the center of the dipole relative to the 
center of the nucleus and the integrand gives the total dipole-nucleus cross 
section for a fixed impact parameter. The nuclear profile function is labeled 
by T A (b) [19]. 



In Fig. 1, the structure functions F^ c (x,Q 2 ) (left panel) and F^ c \x,Q 2 ) 
(right panel) are presented as a function of Bjorken variable x for distinct bo- 
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Fig. 1. The structure functions FJ? C (x,Q 2 ) and F^ c (x,Q 2 ) as a function of x for 
distinct boson virtualities (see text) using the b-CGC model. 

son virtualities using the b-CGC model. It is found that the x-dependence 
is approximately power-like with a effective power which growths on Q 2 , 
\ = A(Q 2 ). Namely, it goes from X(Q 2 = 1 GeV 2 ) ~ 0.12 up to X(Q 2 = 
M|) — 0.224. Similar behavior is verified for the Fi structure function, where 
a flattening is seen at very small-x. This is explained by the QCD evolution 
which is reproduced in the b-CGC model due to a running effective anoma- 
lous dimension. It is interesting to call attention to the unusual behavior of the 
structure function in the limit of large Q 2 and large x, which is more evident 
for Fl- In order to estimate the uncertainty from the theoretical side, we com- 
pute the structure function using the simple GBW saturation model, which 
is presented in Fig. 2. The GBW model does not include the QCD evolution 
in the dipole cross section and in principle it is not suitable for very large 
virtualities. However, the effective power is similar to the b-CGC model with 
F L being distinct at Q 2 = Mf (it could be a consequence of missing QCD 
evolution in dipole cross section). It seems also that the flattening in F L is 
stronger in b-CGC than in GBW. 

In Fig. 3 we present the structure function F 2 as a function of Q 2 for fixed x. In 
left panel one has the virtuality dependence for b-CGC and GBW models for 
two representative values of x (x = 10~ 2 and x = 10~ 6 ). It is verified the usual 
logarithmic dependence as the dipole models reproduce the scaling violation 
for structure functions. There is a small deviation on the overall normaliza- 
tion (a few percents), which is more sizable at large Q 2 . In right panel, one 
has the flavor contribution to structure function. For sake of illustration we 
consider the GBW model and x = 10~ 3 . For light quarks, (d, s) contribution 
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Fig. 2. The structure functions FJ? C (x,Q 2 ) and F^ c (x,Q 2 ) as a function of x for 
distinct boson virtualities (see text) using the GBW model. 

is dominant over u quark, which is evident from the comparison among their 
electroweak couplings. The charm contribution corresponds to 13-14 % of to- 
tal result at this x value. The bottom contribution is very small, suppressed 
by its larger mass. 

An important feature of neutral current neutrino-nucleon interaction within 
the color dipole picture is the geometric scaling property. Such a scaling is 
a basic property of the saturation physics. It means that the total j*p cross 
section at large energies is not a function of the two independent variables 
x and Q, but is rather a function of the single variable r p = Q 2 /Q 2 a _ t {x) 
as shown in Ref. [20,21]. That IS j (Try *p \ X ■ Q 2 ) = a 7 *p(r p ). Similar scaling is 
predicted for lepton- nucleus scattering at high energies [23]. In Ref. [1] we 
have shown that scaling is present in charged current neutrino scattering. It is 
exact for light quarks, whereas is partially violated by charm contribution cs 
and sc due to different quark masses. We can qualitatively show the geometric 
scaling property in the NC neutrino-nucleon scattering by writing down the 
light-cone wavefunctions (for sake of simplicity, we take gv / 9a ~ 1): 
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Fig. 3. The structure function F% N (x,Q 2 ) as a function Q 2 for both models (left 
panel) and its flavor contribution (right panel). 

The equations above are different from the usual expressions for light-cone 
wavefunctions in electron-proton DIS even in the massless case due to the 
distinct weak couplings. Neglecting the quark mass, one can rewrite Eq. (1) 

as, 



a T;{L) (Z° + N-*X)^JdzJ fdff TiiL) {z)K 2 m (y/z(l-z)f) a dt 



where we have fr{ z ) = [(1 — z ) 2 + <2 2 ](1 — z ) z an d fi{ z ) — £ 2 (1 — z) 2 . Here, we 
rescaled the transverse size r to the dimensionless variable f = Qr and suppose 
a scaling behavior for the dipole cross section, a dip = &dip(Qsa,t r )- This last 
assumption is fullfiled by several saturation approaches, as for GBW model 
where a^p = <t [1 — exp(— r 2 Qg at /4)]. Therefore, we have a scaling appearing in 
the NC boson-nucleon cross section, a tot (Z° + N — > X) = a to t( r p = Ql^tlQ 2 )- 
Similar scaling should be present in the nuclear shown in Ref. [1]. 

In Fig. 4 we present the study for nuclear targets. The structure function F^ c 
is presented as a function of x for fixed values of virtualities (Q 2 = 1 GeV 2 and 
Q 2 = 50 GeV 2 ). We select the following nuclei: Carbon, Aluminum, Iron and 
Lead. We stress that for NuSonG experiment we will have an isoscalar target 
similar to A = Fe as for NuTeV. The nuclear shadowing is more intense for 
heavy nuclei, smaller x and virtualities. In particular, for Iron targets the nu- 
clear shadowing is of order 10% at Q 2 = 1 GeV 2 , whereas it reaches to 15-20 
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4. The nuclear ratio R A = F^f /AF^ C as a function of Bjorken variable for 
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fixed virtualities (Q 2 



1 and 50 GeV 2 ) and several nuclear targets. 



% for Lead. Nuclear shadowing is small at large virtualities even at low-x. This 
fact is also verified in charged current interactions as investigated in Ref. [1], 
where the nuclear ratio has been computed for F^ and xF 3 . Some questions 
are in order here. The information on isospin violation is quite relevant for 
electroweak measurements and isospin symmetry is assumed in the NuTeV 
analysis [6]. Violation of isospin symmetry can be one important source of 
NuTeV anomaly. Therefore, better constraint on isospin violation will be cru- 
cial to the interpretation of new neutrino experiments. The NuSOnG experi- 
ment [7] purposes to measure separately AF 2 = y| F£ c (x, Q 2 ) — F^ c (x,Q 2 ) 
in CC W ± exchange and NC 7/Z exchange processes, which can also con- 
strain charge symmetry violation. As NuSOnG (glass target, SiC>2, is very 
nearly isoscalar) could measure F^ on a variety of targets, this will reduce 
the systematics associated with the heavy nuclear target corrections provid- 
ing new information on isospin violation. In particular, for SiO"2 we have 
Z Q = 8, Z Si = 14, m Q = 15.994, m Si = 28.0855. Using A = Z + iV we 

have (N - Z)/A= {A- 2Z)/A, which yields (N - Z)/A 0.000375 for O 

and (N — Z)/A ~ 0.00304 for Si. Therefore, since NuSOnG will provide high 
statistics F^ measurements for a variety of A targets, this will yield an alter- 
nate handle on the isospin violation and will also improve our understanding 
of the associated nuclear corrections. 



As a final study, we compute the associated charm production in NC neutrino- 
nucleus interactions, o"nc( c c) / 0^1? ■ The total NC neutrino-nucleon cross sec- 
tions as a function of the neutrino energy and atomic number are given by the 
integration over available phase space and read as, 
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where Gf is the weak Fermi coupling constant, is the nucleon mass, E v is 
the incident neutrino energy, Q 2 is the square of the four-momentum transfer 
to the nucleon. The variable y = Eh a d/E u is the fractional energy transferred 
to the hadronic vertex with Ehad being the measured hadronic energy, and 
x = Q 2 j2m^E v y is the Bjorken scaling variable (fractional momentum carried 
by the struck quark). The minimum <3mm ^ 0(1) GeV 2 is introduced in order 
to stay in the DIS region. Here, one considers UHE neutrinos, where the 
valence quark contribution stays constant and physics is driven by sea quark 
contributions. Hence, the xF 3 contribution should be negligible and it will be 
disregarded. Similar cross section can be computed for the charm contribution 
by replacing the respective charm structure function. We quote Refs. [22,24] 
for previous analysis on CC and NC neutrino cross sections using the color 
dipole picture. 

Recently, the CHORUS Collaboration [25] has reported an extraction of the 
associated charm production in NC neutrino DIS at E v = 27 GeV. The experi- 
mental determination found a(cc)/a tot = 3.62^^1 (stat)±0.54 (syst) x 1CT 3 at 
the neutrino energy E v = 27 GeV. We quote also previous determinations, as 
the E531 experiment [26] a(cc)/a tot = 1.3±f;i x 10" 3 for energy E v = 22 GeV 
and the NuTeV experiment [27] giving a(cc)/a tot = GAttl x 10~ 3 at E v = 154 
GeV. We have found theoretically a(cc)/a tot = 2.7 x 10" 2 for E v = 27 GeV 
and a(cc)/a tot = 1.35 x 1CT 1 for E v = 154 GeV. Our results are systemat- 
ically above the experimental measurements, which can be explained by the 
extrapolation to very low energies of the color dipole results for charm contri- 
bution. It is well know that the dipole approaches overestimate the structure 
functions for x — > 1. At high energies our predictions are consistent with other 
approaches using small- a; formalisms. 



3 Comments and Conclusions 



As a summary, an analysis of small- a; neutral current (NC) neutrino-nucleus 
DIS is performed within the color dipole formalism. The structure functions 
F% N , F£ N and the nuclear ratio Ra = F% A / AF% N are calculated. In order 
to investigate the theoretical uncertainties we employed two phenomenologi- 
cal parameterizations for the dipole cross section which successfully describe 
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small-x inclusive, diffractive ep DIS data and exclusive production of vector 
mesons. Nuclear shadowing is taking into account through Glauber-Gribov 
formalism. We have found deviations among the models at very small-rr and 
mostly at large virtualities. It is found that small-rr data show geometric scal- 
ing property for the NC boson-hadron cross section as a function of the scaling 
variable r. Although the results presented here are compelling, further inves- 
tigations are requested. In particular, it would be useful compare the present 
calculations to measurements of neutrino-nucleus structure function in smaller 
values of x than the currently measured in the accelerator experiments. This 
is the reason we have focused on the NuSonG purpose, taht will provide high 
statistics F^ c and F^ c measurements for a variety of nuclear targets. We 
also compute the charm content to the total NC neutrino cross section and 
we show that it is consistent with current experimental measurements. So, we 
have a description of small-x structure functions in neutrino-proton interac- 
tions using a color dipole formalism. Despite the limitations of the current 
approach for the low energy neutrinos, our estimates are comparable to NC 
neutrino-nuclei data for charm production. Our calculation is also consistent 
with other approaches in literature in the high energy limit. 
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